SOBOLEV ORTHOGONAL POLYNOMIALS ON A SIMPLEX 



RABiA AKTAS. AND YUAN XU 



Abstract. The Jacobi polynomials on the simplex are orthogonal polynomi- 
als with respect to the weight function W-y(x) = xj 1 ■ ■ ■ x~^ d (l — |i|) 7d + 1 when 
all 7i > — 1 and they are eigenfunctions of a second order partial differential 
operator L-y. The singular cases that some, or all, 71, • ■ ■ , 7d+i are —1 arc 
studied in this paper. Firstly a complete basis of polynomials that are eigen- 
functions of L-y in each singular case is found. Secondly, these polynomials are 
shown to be orthogonal with respect to an inner product which is explicitly 
determined. This inner product involves derivatives of the functions, hence 
the name Sobolev orthogonal polynomials. 



1. Introduction 

The purpose of this paper is to study the limiting case of classical orthogonal 
polynomials on the simplex when the weight function becomes singular. Let T d be 
the d-dimcnsional simplex defined by 

T d := {x G K d : x 1 > 0, . . . , x d > 0, 1 - \x\ > 0}, 

where \x\ := X\ H h Xd- The classical weight function on T d is defined by 

(1.1) W y (x) :=xf ■■■xl d {l- \x\Y d +\ xeT d , 

where 7^ are real numbers, usually assumed to satisfy 7^ > — 1 to ensure the integra- 
bility of W-f on T d . Let Oy = 1 j J Td W~ / (x)dx denote the normalization constant. 

Let n d = M.[x] be the ring of polynomials in d- variables and let II d denote the 
subspace of polynomials of total degree at most n. When all 7* > —1, the bilinear 
form 

{f,9}~-=(h f{x)g(x)W~ / (x)dx 

JT d 

defines an inner product on the space Ii d of polynomials of d-variables. The orthog- 
onal polynomials with respect to this inner product have been studied extensively 
(cf. [7]). In particular, let V d (W^) denote the space of orthogonal polynomials of 
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degree n with respect to (•,•)„. Then, for all P E V, 



(Wy). 



(1.2) 



d 2 p 



X i X j 



dxidxj 



i=l * l<i<j<d 



d 



dP 

dx t 



+ 5Z(7i + l-(l7l+d + l)asi) 



n(n+ \j\ + d)P, 



where 7 := 71 + • • • + 7<2+i- In other words, orthogonal polynomials of degree n 
are eigenfunctions of the differential operator L 7 . 

If some or all 74 equal to —1, the weight function becomes singular and orthogonal 
polynomials with respect to W-y are no longer well defined. In fact, the Jacobi 
polynomials on the simplex defined via the Rodrigue formula can be extended to 
the case of 7, < — 1. In the case of all 7, = —1, they have been used to study 
certain approximation processes in |14i 116] , but this family of polynomials does not 
provide a complete basis. On the other hand, the equation (|1.2|) still makes sense 
if some, or all, ji = — 1. Thus, it is natural to ask two questions: firstly if (|1.2j) 
still has a complete basis of polynomial solutions, and secondly, when the answer to 
the first question is affirmative, if these polynomial solutions are orthogonal with 
respect to an inner product. 

The main results of this paper answer both questions affirmatively. More pre- 
cisely, for each singular case of some or all 7$ = — 1, we have identified a complete 
basis of polynomials that are eigenfunctions of L-y , which are given explicitly, and 
we have found an inner product explicitly, with respect to which these polynomials 
are orthogonal. The inner product turns out to involve derivatives of the func- 
tions, so that the orthogonality is not established in the usual L 2 (W-y) sense, but 
in a Sobolev space. Such orthogonal polynomials are named Sobolev orthogonal 
polynomials. 

The problems are motivated by a recent study [T3J [TB] of analogue problems 
on the unit ball, for which the classical weight function takes the form W^(x) — 
(1 — ||.t|| 2 )''. The orthogonal polynomials with respect to on the ball are eigen- 
functions of a second order differential operator L^. In [18 , a family of orthogonal 
polynomials with respect to an inner product that involves the first partial deriva- 
tives on the unit ball was studied, which turned out to be the eigenfunctions of L_i 
as shown in [13] . 

The Sobolev orthogonal polynomials have been studied extensively when d = 1 
(cf. [9]). The simplex T d becomes, when d = 1, the interval [0,1] and W n is 
the Jacobi weight function x a (\ — x)" , where we have written a = 71 and (3 — 
72. In this case, Sobolev orthogonal polynomials have been studied by several 
authors ([U EJ 21 [11]). More precisely, the Sobolev orthogonality of the Jacobi 
polynomials {Pn~ N '^} n >Q was studied in [3], and the case {Pn N '~ }n>a was 
studied in [4], in both cases N is a positive integer. The particular case of the 
Jacobi polynomials {Pn l ' 1 '}n>o had been previously given in [IT]. In the case of 
d = 2, the simplex becomes a triangle and the weight function is usually denoted 
by W a ,/3 l7 . An observation in [6] shows that the monk basis for W Q)l a l7 is still a 
basis when 7 = — 1 and a, j3 > —1, and they are orthogonal with respect to an inner 
product that involves the first order derivatives. In contrast to one variable, there 
are only a handful papers on Sobolev orthogonal polynomials of several variables 
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[TUl [T2J [131 OH HE], see also [8]. It should be mentioned that the studies in 
[T7l [18] are motivated by problems in numerical solution of Poisson equations [5] 
and in optics. Given this background, it is somewhat surprising that there have 
been so few studies of the Sobolev orthogonal polynomials of several variables. It 
is our hope that the results in this paper will help to kindle more interests on this 
topic. 

The paper is organized as follows. In the next section we recall results on Jacobi 
polynomials and classical orthogonal polynomials on the simplex, where we will 
also present several new properties on the simplex, including how orthogonal bases 
on the faces of T d arise from some of the orthogonal polynomials on T d , which are 
of independent interest. The main results are stated in the third section, where 
we will state the results for d = 2 first to illustrate the results and illuminate how 
the results are obtained, as the results in T d require unavoidably heavy notations. 
Finally, the proofs of the main results, including several lemmas, are given in the 
fourth section. 



2. Orthogonal polynomials on the simplex 

This section contains background results on orthogonal polynomials that are nec- 
essary for latter sections. After a brief subsection on the classical Jacobi polynomi- 
als, orthogonal polynomials on the triangle are described in the second subsection 
and those on the simplex are developed in full generality in the third subsection, 
and orthogonal polynomials on the faces of the simplex are discussed in the fourth 
subsection. 

2.1. Jacobi polynomials. The Jacobi weight function v a> p is defined by 

v a ,p(x) := (l-x) a (l + xf, a,0>-l. 

The Jacobi polynomials are given explicitly by the Rodrigue formula 
r _ i \ n fin 

(2.1) P^\x) = L_J_(1 - x)- a {\ + x)-P— [(1 - xT +a (l + xr+P] 

and they satisfy the orthogonality condition 

(2.2) c a ,p J Pt P) {x)P^ 0) {x)v a ^{x)dx = hfrft8 n , m , 
where 

T(a + f3 + 2) (Q ;3) _ (a + l) n (/3 + l) n (a + /3 + n + l) 

° a ' P ' 2 Q +^+ 1 I> + l)r(/3 + l)' " ' n\(a + p + 2) n (a + j3 + 2n + l) 

and 6 n ,m is the Kronecker delta. The Jacobi polynomial P„ a '^ also satisfies the 
second order differential equation 

(2.3) (1 - x 2 )y" +[/3-a-(a + /3 + 2)x] y' = -n(a + /3 + n + l)y, 

in other words, it is the eigenfunction of the differential operator in the left hand 
side with the eigenvalue — n(a + /3 + n + 1). 

If one of a and j3 is a negative integer, then the weight function v a _p is no longer 
integrable on [—1,1]. Assume a — —I, I £ N, then we have [TS] (4.22.2)] 

(2.4) (^V^ ) (*)=( n t^(^)'ia ) (*). 
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which is well defined for n > I. In the case of a = — 1, defining Pq (x) = 1, then 
„>o is well defined and they still satisfy the differential equation (I2.3[) for 
all n = 0, 1, . . .. They are, however, no longer orthogonal polynomials in the sense 
of (|2.2j) but, as it turns out, orthogonal polynomials with respect to the following 
Sobolcv type inner product O [11] , 

(f,g) := A/(l)s(l) + y (x + l) ti+1 f{x)g'(x)dx 1 A > 0, (3 > -1. 

If both a = (3 = — 1, then the Jacobi polynomials in (|2.ip are well defined for n > 2. 
Furthermore, any linear polynomial will be a solution of the equation (|2 .3[) when 
n = and 1. In fact, the polynomials Pq 1 '~ 1 \x) = 1, p[^ 1 '^ 1 \x) = x + [i and 

Pi-^Hx) = \{x - l)(x + 1)p£#(x), n > 2, 

satisfy the equation (|2.3[) for all n. = 0, 1, 2, Furthermore, these polynomials arc 

orthogonal with respect to the Sobolev type inner product 

(/,<?) := A 1 /(l) 5 (l) + A a /(-l)s(-l) + ^ f'{x)g'{x)dx 

where Ai and A2 are nonnegative numbers, not both zero, and the constant in 
Pi~ 1 '~ 1 \x) is given by ,i = 

Since the simplex T d becomes [0, 1] when d — 1, it is more convenient to con- 
sider the Jacobi polynomials on the interval [0, 1], which we denote by J%*'^ , and 
normalized as 

4°>fi(x) := (1 - x)- a x-e-jL- [(i - x y*« x «+f>] . 

Up to a constant multiple, Jn°'^\x) = cPn a '^'(2x — 1). These polynomials are 
orthogonal on [0, 1] with respect to the weight function 

WaA x ) '■= C 1 -x) a x p , a,j8>-l 
and they satisfy a second order differential equation 

x(l - x)y" + {f3 + l- (a + (3 + 2)x] y' = -n(n + a + /3 + l)y 

2.2. Orthogonal polynomials on the triangle. On the triangle T 2 := {{x,y) g 
M 2 : x, y > 0, 1 — x — y > 0}, the Jacobi polynomials are orthogonal with respect 
to the weight function 

W a ,p„(x,y) := x a y p {l -x- y) 1 , a,/3,7 > -!• 
More precisely, we define the inner product 

(f>9) a ,p n : = c «,/9,7 J f( x ^y)9(x,y)W a ^^(x,y)dxdy, 
where c a ^^ is the normalization constant of W a ^ a given by 

Let V 2 {W a .jj^) denote the space of orthogonal polynomials of degree n with respect 
to this inner product. Then P g V 2 {W a ^ n ) if (P q) a>0n = for all q g IL 2 ^ and 



SOBOLEV ORTHOGONAL POLYNOMIALS ON A SIMPLEX 



5 



dim V^(W a ,/3,j) = (n+ 1). Among many bases for V^(Wq,/3 i7 ), one is given by the 
Rodrigue formula: For < k < n, 

(2.5) P [ ^\x,y) := [W a ^(x, y)]^ dxkQyn _ k [^ +fc /+"- fc (l - x - yV +n ] . 

It should be noted that the elements of this basis are not mutually orthogonal to 
each other. The triangle T 2 is symmetric under the permutation of (x, y, 1 — x — 
y). Parametrizing the triangle differently leads to two more orthogonal bases of 
Vn(W a ,0,-y)- Indeed, if we define 

(2.6) Q n k (x,y):=P^\y,x) and R n k (x, y) := P^ P ' a \l - x- y, y), 

then both {Rl : < k < n} and {Q™ : < k < n} are bases of V%(W a ,p„). Upon 
changing variables from (12.51) . these polynomials are given explicitly by 



SO*, y) = [w a , P , 7 (x, y)]- 1 [z«+«-y+*(i - x - 



and 

Rt{x, y) = {-If [W a , M (x, y)]- 1 dxk{dy ^ Qx)n . k [/^(l ~ * ~ yV +k x a+n ] . 

For the study of the Sobolev orthogonal polynomials, it is necessary to under- 
stand the restriction of the basis element on the boundary of the triangle. For later 
use, we state the following proposition, which is the special case of Lemma 12.101 

Proposition 2.1. For n — 0,1, ... , 

1. The restriction of P^f a) € \%(W a ,p„) on the line x = is J^' P \y). 

2. The restriction of P^ a,l) E V%(W a ,p„) on the line y = is ji 7 ' a) (x). 

3. The restriction of Pfjl^'°^ G ^nO^a^n) on ^ e ^ ne ^ + 2/ = 1 * s Jn^'^iy)- 

There is another basis {V k {a / n) : < k < n} of V^(W a ,p n ) so that p£*/' 7) and 

Vjf^f are biorthogonal. This basis will be the special case of d — 2 of the basis 
(|2.11[) in the following section. 

2.3. Orthogonal polynomials on the simplex. To simplify the notation, we 
set, for x € T d , \x\ := x\ + ■ ■ ■ + Xd and X4+1 := 1 — \x\. Furthermore, for 7$ > —1, 
1 < i < d + 1, we set 

1 ■= (li,---,ld,ld+i) = (l,ld+i) with 7 := (71, . . . ,7^). 
Then the weight function W~ in (|1.1[) can be written as 

(2.7) W^{ X ) = x 1 {1 - \x\)~< d+1 . 

If all 7, > —1, then W-y is integrable and we can consider orthogonal polynomials 
in the space L 2 (W- r ,T d ) with respect to the inner product 

(2.8) (/,S> 7 :=Cy/ /(aOsfaOW^aOda:, 
where c T is the normalization constant of W-y given by 

r(|7l + d + i) 



-v := 1 / / W~f(x)dx — 
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Let V^(W-y) denote the space of orthogonal polynomials of degree n with respect to 
this inner product. As in the case of two variables, there are several distinguished 
bases of V!~(W-y). First we state the basis defined via the Rodrigue formula. 

Lemma 2.2. For n e Nq, 7 G R d+1 and x G R d , define 

(2.9) P2(x) := x-^tl - b|)-T««+i4— M +n (l - M) 7d+1+|n| 

ox n L 

where = ~ n ° „ n d ■ -^f 7i > — 1j 1 ^ * ^ d + 1, £/ien PQ are orthogonal 

1 d 

polynomials with respect to W-y and {P^ : |n| = n] is a basis ofV^W^). 

This is a classical result, see, for example, [Up. 49]. As a basis of V„(H%), these 
polynomials are eigenfunctions of the differential operator L-y in (|1 .2[) ; that is, 

(2.10) L 1 P2 = -n{n+\ 1 \ + d)P2, |n| = n. 

In particular, the differential equation has a complete set of solutions consisting of 
(" + n _1 ) linearly independent polynomials of degree n. 

Remark 2.1. Analytic continuation shows that (|2.10p still holds if some or all 7, 
are < — 1, in which case, however, {P^ ■ |n| = n} no longer contains a complete set 
of solutions. 

If some components of 7 are negative integers, the polynomial P 7 still satisfies 
some orthogonality with respect to a different weight function. More precisely, we 
have the following lemma. 

Lemma 2.3. Let k be an integer, 1 < k < d, and 7 G R d ~ fc+1 with ji > — 1 for 
all i. Let 7 = (7, — m), where m G N fe . Then P^ is orthogonal to polynomials of 
degree at most |n| — |m| — 1 with respect to W^ t o h ; that is, 

f P^-™Xx)Q(x)W J , Oh (x)dx = 0, VQ G nf n| _ |mhl . 

JT d 

Proof. Let us denote temporarily x = (xi,X2) with xi G R d_,c+1 and X2 G K. fe_1 
and, similarly, n = (111,112). Furthermore, let m = (m',mo) with m' G N 
Then, by the Rodrigue formula (|2.9[) . 



J:= / p(?>- m \x)Q(x)W li0k (x)dx 



= J - |x|)™°|^ {x^+^- m '(l - |^|)l»l— o| Q{x)dx 

We perform integration by parts on the last integral. Since 7; > —1, the integra- 
tion by parts over the first d — k + 1 variables carries through without problem. 
For integration by parts in Xd-k+2, the first component of X2, the part that is 
being integrated out, is zero since, if j < mi, then d - _ ; — [x™ (1 — \x\) m "Q(x)] 

® X d—k+12 

vanishes for Xd-k+2 = and \x\ — 1, whereas if j > mi, then for £ = nd-k+2 — j, 
|x" 2 ~ m (1 — |x|)l n l~ mo j vanishes on the boundary. Consequently, we con- 



s' 



elude that 



J = (-l)H J ^ x^+^- m '(l - |x|)l"l- m «|^ [x™'(l - \x\) m °Q( 



dx = 



/ r rd 

since the term insider the bracket is a polynomial of degree |m| + degQ < n - D 
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It should be noted that, instead of 7 = (7,— m), we can assume that 74 is a 
negative integer for i in a subset of{l,2,...,<i+l} in the above lemma. The result 
will be similar and the precise formulation should be clear from the above proof. 

The simplex T d is symmetric under the permutation of (x\, . . . , Xd, Xd+i). This 
symmetry carries over to the Rodrigue basis. Let 

Z d := {l,2,...,d}. 

For a subset S of we denote by \S\ the cardinality of S and by S c the complement 
of S, that is, S c := Z d \ S. 

Let Qk be the permutation group of k elements; its action on a function g : 
R k H> M is denoted by g(xa), a G Qk- Just as in the case of d = 2, sec (|2.6|) . we 
can permute variables and parameters of P^ to obtain different bases of V„(W-y). 

Definition 2.4. Let S = {ii, . . . be a subset of Zd+i and S c = {id+i}- Let 
Is ■= {lS,1i d+ i), where 75 := (7^,. . -,j ld ), and x s = (x n ,.. .,X id ). 
1. If d+l £ S, then S c — {d + 1}, x s = (xi,..., x d )a and 75 = (71, . . . ,7d)cr for 
some a G Q d , and ~/ s = (7s, 7^+1). We denote 

PZ s (*s) = ^ s 7s (i - l^l)" 7d+1 |^ {4 s+n (i - ksir^+w} . 

5. If d+l € S, then S = {d + 1, . . . , id}® for some a G Qd and S c = {id+i}, 
where ij G Z d . Let x s = (x d+1 ,x i2 , . . . ,x id )cr and j s = (j d +i, li 2 , ■ ■ ■ , lijv- 
We define then 

P2 s (x s ) = xZ 7s x; Jid+1 ^- \x1 s+a xj id+1+lnl ) , 

w/iere 9x s , CT = {~d id+l , 9i 2 ,i d+1 , . . . , d id . id+l )o- with di := dxi and 9y := 9, - <9j. 

In the case oi d+l e S, the notation P2 s (xs) agrees with that of (|2.9[) , and it 
is a simple permutation of P£(x). For d = 2, see (|2.6[) . 

Lemma 2.5. Let S and P2 S (xs) be defined in the Definition \2.4\ Then {P2 S (xs) '■ 
n| = n} is a basis of V^(Wy). In particular, P2 s (xs) are eigen] 'unctions of the 
differential operator L-y given by (jl.2p . 

Proof. In the first case, d+l (£ S, changing variable x <— > xa in the integral of 
(Pn,g)^ shows immediately that Pn s (xs) G V^(W-y). 

The second case, d + 1 G S, requires a bit more work. Up to a permutation 
in Qd, we can assume that S = (d + 1, 1, . . . , d — 1). Then S c = {d}, xs = 
(x d+ i,xi, . . . , x d -i) and 7s = (7d+i> 7i: • • • > ld-i)- Making a change of variable 

Ul = X d +1 = 1 - \x\, U2=Xi,...,Ud = x d -i 

so that 1 — \u\ = Xd and, using superscript for the variable with respect to which 
the derivative is taking, 

a l — "ci ' u 2 — "1 ' ■ • ■ ' u d — u d-l u d ' 

we obtain from the Rodrigue formula (|2.9[) that 

P2 s {u) = (-irx^x^ dlnl k s+n *r |n| ] - 

a d ' ' ' U d-X,d 



s 
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Furthermore, under this change of variables, W- i {x) — W~ is {u). Consequently, 
P2 s {x s )g{x)W~ i {x)dx= [ P2 s {u)g(u 2 ,...,u d ,l-\u\)W~ ts {u)du = 

T d JT d 

by the orthogonality of P%, so that P^ s (x s ) € V^(W 7 ). □ 

Another basis of Vn(W-,) is {V2(x) : |n| = n,n G Nq}, where V£ are monomial 
orthogonal polynomials defined by 

'rii + l) n ,(h\ + d)n+|m| m 



(2.11) ™:=SM)" +M n(:')H 

m<n i=l x / v ' 



+ l)r rai (|T + rf )n+|n| ' 



These polynomials satisfy V^(x) = x n + Q n (x), Q n £ U*_ t and VJ G V*(W y ), so 
that V2 is the orthogonal projection of x n in V4(Wy). Furthermore, we have the 
biorthogonal relation 

■ P2( x )vz( x )w^x )d x = nti(y + lUy i + Di-i n! , 

T d (|7| +«+ l)2|n| 

We note that is well defined even if 7^+1 = — 1, since 7^+1 appears only in |— y | 
in the right hand side of (|2.11|) . For latter use, we state one more property of VQ . 
Let ex,...,ea be the standard basis of K d , that is, the ith. coordinate of Bj is 1 if 
i = j, if % ^ j. 

Lemma 2.6. For 1 < i < d, 

(2.12) ^(x)=n i vtT'" d+1+1 \x). 
Proof. We can assume i = 1 and write as 

mi=0m'<n' 1=1 x / w /<»vim />*-r 



l"l 



where x = (x 2l ■ ■ ■ , x^), n' = (712, • • • , «d) and m' = (m2, • • ■ , m<j). Taking deriva- 
tive with respect to xi and shifting the summation index over mi, we obtain 



—vtm = V V ( n"+i m iTT ( ni \ (^ + 1 )"'d'yl + d )"+M 

fel jkiJ?n> fi W (7* + 1)^(171 + d)n+|n| 



mi— 1 m' 



mi —1 m' <n 

"e E (-ir + 'H-vf ni - 1N ) -fa±ik 

„&J<V V mi ;(7i + l) mi+ i 

x yr fni\ {ji + 1 )n,(l7l +rf)n+|m| + l x ,, 



i=2 V 71 *' + 1 )"»id7l + d )n+\n\ 

Since (a) n = a(a + l) n _i, it holds that 

(Ti+jjni ((7i + l) + l) ni _i (|7l+d) n +|m|+i _ (l7l+rf + 2) n+ | m |-i 

(7l + l)mi + l ((7l + l) + l)mi ' (M + tOn+H Cl"rl + d + 2 W|n|-2 ' 

from which (|2.12p follows. □ 
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The relation (|2.12|) allows us to take higher order derivatives of VJ. One inter- 
esting consequence of the lemma is that are orthogonal polynomials for a family 
of other inner products that involve differentiation. To state the result, we need 
further notation. For a subset S = . . . , ij} of Z<j, we define 

dx s dx{ 1 . . . dxi j 
Lemma 2.7. Let ji > — 1 for 1 < i < d + 1. For m < d, define 



(2.13) 



|S|=j 

where As are nonnegative numbers. Then the polynomials V£ in (|2. 1 1[) are orthog- 
onal with respect to the inner product [•, •]-,,. 

Proof. By (|2.12[) . taking repeated derivatives of V^, we have 

with S = {i\, . . . , ij}, which implies, by the orthogonality of V^J^,. 1 + J ,7<i+1+ ^ ( 

»x ••• la- 
titat is orthogonal with respect to 

and, consequently, orthogonal with respect to [ m ,-]y Since all As > 0, [■, -]"y is an 
inner product. □ 



Remark 2.2. Using the relation (|2.12[) repeatedly, we can also include partial deriva- 

dx. 1 ■■■dx 



tives of higher orders, such as m f m fc for m = (mi, . . . , mfe) with fc < d, in the 



definition of [f,g]y For our purpose, the definition in (|2.13l) is sufficient. 

A couple of remarks are in order. We note that [•, •]§ is not an inner product by 
itself since, for example, [f,f]s = for f(x) — 1. Even though {VJ : |n| = n} is 
a basis of V„(W T ), not every basis of V„(W-y) is orthogonal with respect to [-, 
For example, when As are not all zero, P^ is not an orthogonal polynomial with 
respect to [•, 

2.4. Orthogonal polynomials on the faces of the simplex. In addition to 
the standard basis ei, . . . , e<j of R d , we further set eo = (0, . . . , 0). The points 
eo, ei, . . . , e<j are the vertices of the simplex T d . The boundary of the simplex T d 
can be decomposed into lower dimensional faces. 

Definition 2.8. Let S be a subset of Z^+i with 1 < \S\ < d. Then 

Ts HS| := {x £ T d : Xl = 0, % e S} 

is a (d — \S\)- dimensional face of T d . 

For 1 < \S\ < d, the number of (d — |5|)-dimensional faces is ( d \s\)- In particular, 
for \S\ = d, the 0-dimensional faces are precisely the vertices eo, ei, . . . , ej of T d , 
whereas for |5| = 1, the (d— l)-dimensional faces are , j = 1, . . . , d, d+1. Each 
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face Tg S is a (d — ISD-dimensional simplex and, moreover, Tg S is isomorphic 
to the simplex T d ~l s in variables {xi : i 6 S c }. Furthermore, the restriction of Wy 
on Tg 5 is a weight function of the same type in variables {x{ : i € S c }. 

Definition 2.9. Let S = {i\,...,id} be a subset of Zd+i and S c = {id+i}. Let 
Is '■= {ls,7i d+1 ), where 73 := (7^, . . . ,Ji d ), and x s = (x tl , . . . , x id ). Let P2 s (x s ) 
be defined as Definition \2.4\ 

1. Let d + 1 ^ S. For Sj C Id with \Sj\ = j, 1 < j < d — 1. define 

H ts 3 ( w i) ■= span{PJ(a;) : |n| = n, and n e = 0, W € Sj}. 

2. Let d + 1 e 5. JTiera S = {d + 1,4,.. . twi/i {4, • • ■ Jd-i} C Z d and 

\ S = For = {d + 1,4, • • .,4-i} with 1 < j < d-1, where 

Si = {d + 1}. define 

Ks,( W -y) ■■= spMP^'^^s^xs^) :\n\=n with n e = 0, 1 < £ < j}. 
We further define H* >s<t (W y ) = and H^ Si+i (W^) = 0. 

Lemma 2.10. The space iJ„ s-^Wi) satisfies the following properties: 
00 H n,s 3 ( w ~/) is a subspace ofV*(W y ); 

(ii) The restriction of g.(W-y) on the face Tg~ 3 is the subspace V^ _J (W-y| T <i-j), 
where W-y\ T d-j denotes the restriction ofW-f on the face Tg~ J , that is, 

Ks 3 ( Wn ) I T *-i = Vt 3 (W n I T « ;j ) 
in the variables {xi : i G Sj}. 

Proof. The item (i) follows immediately from Lemma 12.51 For the proof of item 
(ii), we need to consider two cases, d + 1 ^ S and d + 1 € S. 

Case 1. d + 1 \t S. Then S = Let Sj = . . . ,ij} C Z^, Sj = {ij+i, • • • , id}, 
7s? = (Tij+u- • • >7iJ and n S c = (m j+1 ,. . -,n id ). U n e = for I £ Sj, then in 
()2.9|) can be written as 

_ ,Jn s c| 

^(x)=^ ? (l-H)- w ^i 
which, when restricted to the face Tg~ 3 , becomes 

_ „|n s e| 
P2(x)\ Td - D = Xg^Hl ~ IsecD-W+l—L 

dXgJ 

precisely the Rodrigue formula in the variable xs c for the weight function W~| a-j 

and, as such, is an element of V„~ 3 \Wy\ T d-j) ■ Furthermore, the collection of such 

polynomials forms a basis of the latter space. This completes the proof of (ii) in 
this case. 



7s c + n s? 



3 (i-M) 



7d + i + |n s c| 



7 s c+n S c 



1d + l-\-\n S c 
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Case 2. d+1 E S. Up to a permutation in Qd, we can assume S — (d+1, 1, . . . , d— 1). 
Then S c = {d}, 7 S = (~fs,7d) with 7s = (jd+i, 71, . • . , Jd-i) and Sj = {d + 
1,1, . . . ,j — 1}. Making a change of variable 

Ui = X d+1 = 1 - U 2 = X X ,... ,U d = Ed-l, 

we have 1 — |tt| = ce<j, W T (a;) = W Ts (u) and Pn s (^s) = P2 s {u), as in the proof of 
Lemma [2751 so that 

Hi iS .(W y (x)) = Hl s ,(W^ s (u)) 
where S 1 * := {1, 2,...,j} and 

< s; (W.y») = span{P^( W ) : |n| = n, n e = 0, W E 5*}. 
From the case that d + 1 £ S, it follows that 

<s; (W 7s )| = V*"' (V 7s I jid-j ) , 

which is equivalent to 

< s 3 (W y )\ Tl -i = (W, I T , ;3 ) • 
This completes the proof. □ 

Lemma 2.11. Let H% s (W~f) be defined as in Definition \2.(A Then H^ S ,{W^) is 
a subspace of V„(Wy*), where 7* = (7* , . . . , * s defined as follows: 7^ E R 

are arbitrary for I E Sj and 7| = 7j otherwise. In particular, every element of 
s .(W~f) satisfies the equation L^.P = — n(n+ (7* | + d)P 7 where 7^ are arbitrary 
for'i&Sj. 

Proof. In the Rodrigue formula PJ of (|2.9p . if = 0, then the derivative with 
respect to xi does not appear, so that, after canceling out xj ll! and xj e , does 
not appear in which means that we can assign 7^ any value. In other words, 
P2(x) = 0) ifn e = 0,£€ Sj. Consequently, H^ s . (W 7 ) E V*(W*). □ 

3. Main results 

Our main results on the simplex are fairly involved in notations. We shall state 
our results first in the case of triangle, that is, d — 2, and use it to explain how we 
arrive at the results. 

3.1. Sobolev orthogonal polynomials on the triangle. As stated in the intro- 
duction, we are facing two problems. The first one is to find a complete solutions 
of polynomials for the differential equation (|1.2p in which d = 2 and 7 = (a, /3, 7). 
To facilitate the discussion, we define the following subspaces of V^(W ail a 7 ), 

H n ,i(w 7tl3 ) := span{P (0 / ,7) (a;,2/)}, H na {w lta ) := span{P (0 „"' 7) (y, x)}, 

H n> 3(w at p) := span{P (0 /' Q) (l -x- y,y)}, 

which can be formulated as special cases of the spaces in Definition 12.91 According 
to Proposition [27Q H n ^ 1 (w 1 ^)\ x=0 = V*(w^p), H n ^{w 1 . a )\ y=0 = V^(w^ <a ) and 
Hn,a{ w a,p)\x+y=i — V\{w a ,fi)- In the following we adopt the convention that 

V^(W Q)/ 3, 7 ) = and H n . 3 = 0, if n < 0. 
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Theorem 3.1. For n € No and at least one of a, ft, 7 being —1, the differential 
equation 

L a ,p,-< u = —n(n + a + ft + 7 + 2)u 
has a solution space U^{W a ^ t j) of polynomials of degree n, which has dimension 
n + 1 and can be decomposed as a direct sum as follows: 

Z^(W qA _i) := (1 - x - 2/)V^_ 1 (W a , (3 ,i) + ff n , 3 K ? ), n > 0, 

Z£(W a ,_i,_i) :=y(l -x- y)V 2 n _ 2 (W aX1 ) + yH n - 1)3 (w ail ) 

+ (1 - x - j/)ff n _i )2 (u)i ia ) 

for n > 1 and Wq (W ai _i._i) := span{l}, and 

^(WLi,-!,.!) :=xy(l - x - |/)V^_3(Wi, u ) + a V ff B _ 2 , 3 (t« 1 ,i) 

+ i(l - a: - y)i?„-2,2(wi,i) +2/(1 - a: - y)H n _ 2A (w 1 . 1 ) 

for n > 2, Wq(M / _i, -1,-1) := span{l} and U\ (W-\^\-\) :— span{x + C\, y + c 2 }, 
where C\ and c 2 are arbitrary constants. 

From the decomposition, we can deduce further relations among the subspaces. 

Corollary 3.2. The subspaces in the Theorem \3.1\ satisfy the relations 

Z^(W a ,_i,_i) =yUl_ 1 {W a , lt _ 1 ) + {l-x-y)H n _ 1 , 2 {w lta ), n>\ 

Z^(W_i,_i,_i) = 3^-1(^1,-1,-1) +y(l- x - y)// n _ 2 ,iKi), n > 2. 

Remark 3.1. In the statement of the theorem, we could also assume either a = — 1 
or ft = — 1 instead of 7 = — 1. The result then can be obtained via permutation 
of (a, ft,"]) and the corresponding permutation of (x, y, 1 — x — y). We give one 
example explicitly: 

W^(W_i,_ 1)7 ) = xyV r 2 l _ 2 (W / i,i,7) + xHn-x^w^s) + j/JJ n _i ) i(w 7) i). 

Our second problem is to identify an inner product under which the eigenspaces 
in the Theorem 13.11 are orthogonal. 

Definition 3.3. Let Xi and A.y be nonnegative numbers. Define the bilinear forms 

(/' 9) a 0-1 '■= / [xd x f(x, y)d x g(x, y) + yd y f{x, y)d y g(x, y)] x a y l3 dxdy 
Jt 2 

+ Ai [ f(x,l- x)g(x, 1 - x)x a (1 - cc^da; 
Jo 

(/,.9) Q -1 -1 : = / d y f(x,y)d y g(x,y)x a dxdy, 
Jt 2 

+ Ai / 5 !C /( : E,0)9 x 5(^0> a+1 ^ + Ai,o/(l,0)ff(l,0) 
Jo 

(f,9)-i -1 -1 : = / d xy f{x,y)d xy g{x,y){l-x-y)dxdy, 

JT 2 

+ Ai / ^/(x, 0)^5(21, 0)d.x + A 2 / d y f(0,y)d y g(0,y)dy 
Jo Jo 

+ Ai, /(1, 0)5(1, 0) + A ,!/(0, 1) 5 (0, 1) + A ,o/(0, 0). 9 (0, 0). 
Theorem 3.4. Lei a, ft > — 1. For n = 0, 1, 2, . . i/ie following holds: 
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(1) If Ai > then (f,g) a a _ x is an inner product for which lf%(W a ,p,-i) is the 
space of orthogonal polynomials of degree n; 

(2) If Ai > and Xi q > 0, then (/, g) 1 _ 1 is an inner product for which 
W„(W a _i_i) is the space of orthogonal polynomials of degree n; 

(3) IfXi > 0, A2 > and at least one of Ai,q, Ao,i, Ao,o is positive, then (f ', g) _ 1 _ 1 _ 1 
is an inner product for whichli£(W—i —i —i) is the space of orthogonal polyno- 
mials of degree n, where we redefine L/%(W—i t —i —i) as 

^1 (W_i,_i,_i) = span <^ x - - — — — ,y 



Ai,o + Ao,i + Aq,o Ai,o + Ao,i + Ao,o 

With respect to (•, ■) a a _ v the decomposition U%(W a ,p,-i) is an orthogonal one, 
that is, 

U^w a , p ,-i) = (l-i-i/^^AijQ^K?), 

as will be shown in (iii) of Theorem 13.71 

Again, by symmetry, we can state the inner product if other parameters are — 1 . 
We give one example. 

:= J J,dy ~ d x )f(x,y)(d y - d x )g(x,y)(l -x- yfdxdy 

+ Ai / d x f(x,0)d x g(x,0){l-xy +1 dx 
Jo 

+ A 2 f d y f(0,y)d y g(0,y)(l - yV +1 dy + A , o /(0, 0) ff (0, 0) 
Jo 

where Ai, A2 > with at least one of them positive and Aq^ > 0. 

Remark 3.2. Notice that {f,g)_ 1 _ 1 contains two terms of one-dimenional inte- 
gral, whereas {f,g) a _ 1 _ 1 contains only one such term. Evidently, by symmetry 
we could include another term in {f,g) a _ 1 _ 1 - In fact, it will become clear in the 
discussion below that we could include even more terms if we are willing to use 
higher order derivatives. Apart from consideration of symmetry, as in the three 
vertices in (/, <?}_ 1 _ 1 _ 1 in Definition [3T31 we shall keep the number of terms in the 
inner product minimal below. 

These theorems will be special cases of the results on T d . We will not give 
separate proof for d — 2. Instead we give an indication below on how the main 
terms of the inner products are identitified. We start from an observation in [6] that 
the monic orthogonal polynomials vj. a ^'' 1 \x, y) are also orthogonal with respect 
to the inner product, when a, /3, 7 > — 1, 

(3-1) [/, g] a ^ n = (/, g) a ^ n + Ai {d x f, d x g) a+i pn+1 + X 2 (d y f, d y g) a p+l l+1 , 

which is a special case of Lemma [2.71 Since V^f' 1 are well-defined if 7 = — 1, 
letting 7 1— > — 1 preserves the orthogonality, which is how the orthogonality in 
(f,g) a p _j was established [B], but the method does not give the decomposition 
of the orthogonal space. Indeed, the second and the third terms in the right hand 
side of [f,g]a,/3,j make sense for 7 = —1 as well, whereas for the first term we 
parametrize the integral over T 2 as 

fl rl-y 

f(x,y)dxdy= I / f(x,y)dxdy, 



t 2 Jo Jo 
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and apply the limit relation 

lim — |- / f(x)x»dx = /(0) 

M ->-l+ 1-1+ I Jo 

to conclude that 

lim [f,g] a ,p,*f =c a ,p \ f(x,l-x)g(x,l-x)x a (l-x) l3 dx 

7->- 1+ Jo 

+ Xi(d x f, d x g) a+10 o + \n(dvf,d v g) aiP+li0 , 

which is a constant multiple of (f,g) a a _ 1 when A2 = Ai. Since vjf^' (x) are 
well defined for all < k < n, their orthogonality is preserved under the above 
limit. 

The above limit process can be continued one more time by taking, for example, 
(3 — > — 1, which gives 

lim J/,sk,3, 7 = c{f,g) a _ 1 _ v 

where c is a constant. However, since VL"' ^ (x) are not well defined when k = 0, 
the limit process at this stage no longer gives the complete answer for our problem. 
If we take the limit a —> — 1 in {/, g) a _ l _ 1 , we end up with 

lim (/, g) a _ x _! =Ai / ^/(x, 0)d x g(x, 0)dx+ 
<*->- 1+ ' ' Jo 

+ A 2 /" fl v /(0, j/)^ 5 (0, + A li0 /(1, 0)<?(1, 0) 
Jo 

which however is no longer an inner product since, for example, it is equal to zero 
if f{x, y) = g(x, y) = xy. In order to obtain the main terms of the inner product in 
(fid)-i -1 -1 we need to go one level up by noticing that, by Lemma V^"' /3,7 ' 1 
are orthogonal with respect to yet one more inner product: For a, /?, 7 > — 1, 

(3-2) [/, g] a , 0a = (/, g) a , Pn + Ai {d x f, d x g) a+l f3 l+1 

+ h(d y f, d y g) aj3+l l+1 + \i,i(d xv f, d xy g) a+1 :(3+ i i7+2 

where Ai,A2,Ai j i are nonnegative numbers. Taking the limit three times in the 
above inner product leads to the main terms of the inner product (f,g}_ 1 _ 1 _ 1 - 

We emphasis, however, that the above limiting process only suggests possible 
form of the Sobolev inner product. It neither gives the final answer nor it proves 
the orthogonality. 

3.2. Sobolev orthogonal polynomials on the simplex. For a fixed positive 
integer k, 1 < k < d, we consider with = ... = 7d+i = — 1. For 

convenience, we introduce the notation 

where lk ■= (1, . . . , 1) £ Z fc and — lk = (— 1, ■ ■ ■ , —1), and we always assume that 
the subindex of W-y is a vector in R d+1 so that, with the length of the vector — lk 
being k, the length of 7 is d — k + 1, 7 = (71, ... , jd-k+i)- In particular, 

ww*) = xf . ..xf_- k k n^- k+ 2 ■ -^(i - Mr 1 - 
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In these notations we can allow k = d + 1 if we regard 7 £ M. d ~ k+1 as nonentity 
when k = d + 1 and write W-i d+1 (x) = W-i(x). 

Recall the subspace H* s ,(W*y) in the Definition 12.91 From Remark 12.111 since 
the values of 7$ for i £ Sj can be arbitraty in the basis of H% g.(W^), we can set 
{7i = 0, i £ Sj}. Let 

Sk,j(<r) ■= ■ a = 0, (ei,. ■ • ,e d+1 ) = (7, (0 fc _ i; 

so that Sk,j(cr) contains the indices of zero elements in (Ofc_j, lj)o~, and we treat 
Skj(cr) as a vector whose elements are arranged according to their order in (7, (Ojk-j, 
To simplify the notation, we define for j = 0, 1, . . . , k — 1, 

where if a = id, then S/^o = {d + 2 — fc, . . . , d + 1} and lo is a nonentity. It then 
follows by Lemma \2. 101 that 

(i) H*(W^ 0h _ jilj)a ) is a subspace of V^{W l! ( 0k _ j>1 - )rT ); 

(ii) The restriction of H^(W 7 ^o kjlj ^ rT ) on the face T g ~^~^) agrees with the 
space V, < ^ (fc ~- j) (W 7! ( 0fc _.. l3 ) (T | T d-( fc - 3 )) with variables in {cc; : i £ (S k ,j{cr)) c }. 

In particular, if a = id, then Sk,j = {d — k + 2, . . . , d + 1 — j} and 

^n(W / 7,(0 fc - J -,l J ))| T d -( fc -J) = V^~ (fc_j) [W( 7: o k _. A .)\ T d-(k- 3 ) J 

with variables xi, . . . , Xd-k+i, Xd+2-j, ■ ■ ■ ,x d . Furthermore, by definition, we also 
have (W 7 ,oJ = if fc = d or d + 1. 

We adopt the convention that X^=j = if fc < j regardless the values of a,. 

Theorem 3.5. Let k be a fixed integer, 1 < fc < d + 1 and 7 G For 

neNo, the differential equation 

L-y —i k u = —n(n + I7I + d — k)u 
has a solution space Ll*(Wy —t h ) of polynomials of degree n with 

dim^(W 7 ,_ 1 J=( n + ^" 1 ) 5 

which can be decomposed as a direct sum as follows: 

(i) if k = 1, then Uq(W 1 -i) — span{l} and for n > 1, 

(3.3) Un(W 7 ,-i) = x d+1 V*_i(W^i) + H*{W lfl ). 

(ii) if 2 < k < d, then Uq (W-y—i h ) — span{l} and for n>l, 

(3.4) U*(W lt - lh ) =x d+1 - {k _ 1 yx d+1 Vt k (W J , lk ) 

k-l 

+ E E x jA<r)BLj (W^o^y) +H*(W^ 0k ), 
where x d+ i = 1 - and with x := 1 and (ei, . . . ,£ d+ i) = (0 d - k +i, (Ok-j, lj)o), 

d+l 

x-j,k{o~) '■= Y\_ x $i with Si = i if £i = 1, and Si — if Ei = 0. 

2=1 
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(Hi) if k = d+1, then Uq(W-i) = span{l} ; Uf(W-i) = spanjxi +ci, . . . ,Xd + Cd} 
for arbitrary real numbers a,. . . ,Cd, and for n>2, 

d 

U d n {W_i) = Xl ■ ■■x d+1 Vf l _ d _ 1 {W 1 ) + E x i4+^)HLj (W {0d+1 _ jilj)<r ) , 
where, with (e u . . . ,e d +i) = (O d+ i-j,lj)a, 

d+l 

x j.d+i( cr ) '■= Y\_ Xs i &i — i if Si — 1, and 5^=0 if = 0. 

i=l 

In the case of k = d + 1, that is, all -fj = — 1, the differential equation becomes 

L_iP = -n(n - 1)P. 

For n = 1, any polynomial of the form P,(a;) = a;, + Cj with Cj being an arbitrary 
constant satishcs the equation L_ 1 P = 0. The values of the constant will be fixed 
when we discuss the Sobolev orthogonality. 

As an example, we give the decomposition for k = 3 and k < d explicitly: 

Un( W l,-h-h-l) = Xd-lXdXd+lV^iWj.i,!.!) +x d -ix d H%_ 2 (W 7 ,l,l, ) 

+ x d _ix d+l H*_ 2 (W 7) i i0 ,i) + x d x d+1 H*_ 2 (W 7i0 ,i,i) 
+ x d -iK-x (W 7 ,i, 0l o) + Xdttt-x (W 7 ,o,i,o) 
+ x d+1 H*_i (W M1 ) + H d {W M0 ). 

Our next theorem shows that the elements in U d (W lt -i k ) are orthogonal polyno- 
mials with respect to an inner product of the Sobolev type. We need the following 
notations: for 1 < k < d — 1, define 

x fe := (x!,...,x k ), |x fe | = xx + ... + x k , 7 fe := (7i,...,7fe). 

Recall that, for I C Z d +i with |7| = i < d, the face if - * of the simplex T d 
is defined by T/~ 4 = {x G T d : ^ = G /}. If 7 = we define 

x/ . , • • • , Xj i ) . 

Definition 3.6. For 7 G K d , 7^+1 = —1 and A > 0, define 

(f,9)^ - E J Td ^tf dX + X L-, fWtottV - |x,-i|)-^_, 
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For 2 < k < d, let m fe = {d - k + 2, . . . , d} and = {d - k + 2, . . . , d + 1}. For 

7 e R d+1 ~ k and A, A,, A/ > 0, define 

l Td dx m " dx m " 

fe-2 



i=l JCm fc 



TCrrifc 

|J|=i 



i=l 

+ A 



e ^ S T ^ +1 x ^M^ +i - kd " d - k+i 

/ /(xd-fc+i,0) 5( x d-fc+i J °) x Lfc+i dx d-fc' 



where if k — d, we replace the last term by A/(ei)g(ei) ; with A > 0. Finally, when 
all 7j — = — 1, for Xj > and A^.o > 0, define 

f r) d f r) d n 



d-1 



i=l /C 
|J|=i 

+ A ,o/(eo).g(eo) + Ai, /(ei)g(ei) + . . . + X dfi f(e d )g(e d ). 

In the following we shall write (/, <?} _j with k = 1, 2, . . . , d + 1 for all three 
cases. Whenever (/, <?) _ lfc is an inner product, we can define a space of orthogonal 
polynomials. Our next result says that the space of orthogonal polynomials with 
respect to (f,g) _ ± is exactly W^(W 7i _i fc ). First, however, we need to specify 
Uf{W-x): We redefine 

Uf{W-x) = span{xi + ci, . . . ,x d + c d } with Cj = — ^ 3 '° — - — . 

^0,0 + Al,0 + • . • + Ad,0 

Theorem 3.7. F or 1 < k < d+1, let 7 e M d+1 -' i; with jj > -1, 1 < j < d+ 1 - k. 
(i) (/, g) _j is an inner product if A, A;, A/ > for 1 < k < d and, in addition, 

at least one of Aj,o is positive for k = d + 1. 
(ii) i/„(W 7) _i fc ) is the space of orthogonal polynomials with respect to (/, g) _ t . 
{Hi) In the case of k = 1, the decomposition of (|3.3p is orthogonal, that is, 



Remark 3.3. For k > 2, the direct sums in the decomposition of li™(W~ is 
in general not orthogonal sums. In other word, the elements belong to different 
parts of the decomposition (I3.4[) may not be mutually orthogonal in general. On 
the other hand, the Gram-Schmid process is applicable if an orthonormal basis is 
desired. 

It should be mentioned that, as stated in Remark [321 it is possible to add more 
terms or higher order derivatives in the definition of (/, g) _j but we strike to 
keep the number minimal. 
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4. Lemmas and Proof of Theorems 

In the first subsection we prove several lemmas that will be needed for the proof 
of Thcorem l3.5l The proofs of the main results are given in the subsequent sections. 



4.1. Lemmas on Rodrigue basis. Our first two lemmas are analogues of (|2.4j) 
of one variable with one of the element in 7 being — 1. 



Lemma 4.1. Let 7 6 M. d and 7<j+i = — 1- Then, for m > 1, 1 < i < d, 

d 

P^- 1 \x) = (l-\x\)J2- 
Proof. Applying the product rule 



"«(7« +n i ) i3( ^,i) 



— ) [(1 - x)g{x)] = (1 - x)—g{x) - n-^^g{x) 



dx r ' 



dx r - 



multiple times within the Rodrigue formula, we obtain 



#|n| 



P^\ x )= x -i— x^{l-\x\)\- 



in which the first term in the right hand can be written as 



, flN r 



c 7+n (l 



("1 



l+_M x -7^ 

|n| ^ dx n 

l—l 



x 7+n (l- |a;|) 



c^ +n (l - |x|)l n l 



|n|-l 



9 



(xT+"(l-\x\)W) 



so that, putting the two identities together, we conclude that 



0M- 



1 1 i= 1 



c 7+n-e 4 ( 1 _| a .|yn| 
'(*), 



which completes the proof. 

Lemma 4.2. Assume 7, = —1 /or some i, 1 < i < d. TTien for n% > 1, 



□ 



P^'^\x) - -(|n| + 7(i+ i)^™' 7£i+i; W- 

Proof. Without loss of generality, we can assume i = d. Setting 7 
and x = (x.d-i, x d), and applying the product rule 



,(7+2ei,7d + i) , 



(7 d -i,-l) 



d n 

l x 9(x)} = x -^9{x) + n-j-^gix), 



dx r ' 



multiple time with the Rodrigue formula, we obtain 

Q\n\ 



dx 11 



n d x d 2r^-\x\r~< d+i 



x 

dx n ~ e " 



-n i-ed^ _ l^iynl+74+i 

a; 7+n (l - |a;|)l n l+7<*+i 
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Splitting the partial derivatives in the first term as = g x „-<, d q~ ; a quick 
computations establishes the stated identity. □ 

By repeatedly applying the above lemmas, we can deduce a relation when more 
than one ji equal to —1. We state the results in the following corollaries. 

Corollary 4.3. Let 7 e R d+1 - fc . For n d _ fe+2 > l,...,n d > 1, 

p^- lk -^ d+l) (x) H-^W-k ■■■x d 

x JJ(|n| + ld+1 j + l)P^lt^ d i ; l) (x). 

In particular, it follows that 

Xd+2-k ■ ■ ■ 

Combining Lemmas 14. f I and 14.21 , we can also deduce the following lemma: 

Corollary 4.4. Let 7 e R d+1 - k . For n d - k+2 > 1, . . . , n d > 1, 

d 

(4.1) P^-^(x) = x m ■ ■ ■ x d (l - \x\) ]T K^-lf-io,!^)' 

i=l 

where, for 1 < i < d, 

An il7 i = (-1) rn II(N-J)- 

1 1 3=0 

For later use, we need to reverse the expression in (14. ip , which is the content of 
the following lemma. 

Lemma 4.5. Let j — (j'2, . . . ,j d ). Then there exist constant fij such that 

(4.2) x m ■ ..x d (l - \x\)P^(x) = £ Mj^( ( r„Hil ) + i,j) + (o,i fc _ 1 )( :E )' 

j'i<n« 

where the sum is over ji for i = 2, 3, . . . , d. In particular, for n > k, 

x d+2 -k ■ • • x d (l - \x\)V*_ k (W yAk ) C V*(W 7 ,_iJ. 
Proof. To simplify the notation in the proof, let 

X k := x d _ fc+2 • • • x d (l - |x|), n fe := (m, . . . , n k ) 6 N*. 
From (14. ip it follows that 

p(7,-l fc ) — Y, \ p(7,lfc) („A 
r (n,O d - k ,\) ~ yv -fc A n>7l- r ( n -l,0)^''' 

which proves (|4.2I) for n = (n, 0, . . . , 0). Notice that the subindex n of P2 is always 
in Nq ; so that we have suppressed the subindex of 1 = lfe_i for the polynomial in 
the left hand side and the subindex of — Q d -\ for the polynomial in the right 
hand side. We shall keep this convention below. Again by (|4.1j) . we have 

\ p(7,!fc) (~\ 1 \ p(7,lf») f~.\ 

A «i,7i- r (ni-l,n 3 ,0) W T *'»2>7a- t (ni,n2-l,0) WJ ' 

which shows by induction on n 2 that 



^ > (ni,n 2 :O t! _ (fe+1) ,l) 



n 2 

x fe p, (iili, n ,(i) = y w p, (7 'T lfc) ... . n 



J'=0 
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where the fij can be determined inductively, which proves (|4.2[) for n = (m, ri2, ■ ■ ■ , 0). 
Furthermore, the above proof shows, in fact, that 

"2 

^- P ((ni,n2,0 <i _( fc+1) )a,0)( a: ) = ^2 ^ P «ni+n 2 -] + l ,j,0 d _ (fc + 1) )<x.l) ( X ) 

3=0 

for every a £ Ga-k+i, the permutation group of (d — k + 1) elements. In the next 



step, we deduce from (14.1[) that 

jP (ni,n 2 ^Lo £i _ (t;+2) ,l)( a; ) = [^"im-ffni-l.na.ns.O)^) 

+ ^"2.72'f(^'i,n]-l,n3,0)( :r ) + ^s.'W-^ni.na.ns-l.O)^)] - 

For n.s = 1, we can use induction on no to deduce, from XtP,7' „ . nifz), 

the desired (gj) for X h P[^ m (x) and, hence, for XtP^^^^*)- 
Then, by induction on 7J3, we can deduce the desired (|4.2I) for n = (m, ri2, 713, 0). 
It is now clear how to proceed, by induction if necessary, to conclude that (I4.2[) 
holds for n = (n d -k, 0). 

We can go beyond d — k by using (|4.1I) and the fact that A m , T = if m = 1 and 
r = — 1. Indeed, (|4.1I) implies that 



ci-fc+l 

^ ( " Y,_1 * ) ^(as)=-y* V V y .P^' lk) n , (x), 

(n d _ fc+ i,l)V / ™»>7« (n d _fc + i,0)-ei V /> 

j=l 

which allows us to use the induction on n^-fc+i to prove (|4.2I) for n = (n^-fc+i, 0). 
Furthermore, the process can obviously be continued, by induction if necessary, 
until we reach n = (ni, . . . , n d ). This completes the proof. □ 

4.2. Proof of Theorem HU By (03]) in the TheoremEU the space W*(W 7 ,_iJ 
consists of three main terms. Since, by Lemma 14.51 the elements in the space 
Xd+2-k • ■ ■ Xd+{V^_ k (W lt x k ) are orthogonal polynomials of degree n in V^(W 7! _i fc ), 
by Remark [2H1 they satisfy the equation 

(4.3) L 7i ^i k u = —n(n + 7I + d — k)u, 

which takes care of the first term. The second term of the decomposition in (|3.4|) is a 
large sum. Since the permutation is acting on (0fc_j , lj) and the differential opera- 
tor L J ^o k _ j ^ j is clearly invariant under such permutation, we only need to consider, 
for 1 < j < k - 1, the cases of x jtk H^_ j {W lf0k _ ji x j ) with x j>h = x d+2 -j ■ --x d+1 
and Xj.kHf l _j(W~ fy o k j l .i j .o) with Xj t k — Xd+i-j • ■ • x d . We consider the case of 
H^jiW^o^,^) first. Recall that 

//;; f (w y , 0l ., i = ffl iA , 3 .(ww^) with s kJ = \<i /,• • ■>. . . .../ j • i). 

From Lemma 12.111 since the parameters in {7^ : i £ Sk,j} do not appear in the 
basis of -ff^_ J (W 7! o fc _j,i J ), the values of these parameters can be arbitrary. As 
result, we can write 

(4-4) ^t-jCWW,.!,) - ^-i(Wr,-i»-,.ii)- 

On the other hand, from Lemma 14.51 for have 

X d+2 .j---X d (l - IxDVtjiWy,-!^,!,-) C V'fl^.-lJ. 



a 
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Since, by Lemma £321 H^_ j {W 1 -i k _ jAj ) C V^^.-i^,^), together with 
(l4~4l) it follows that 



Xd+2-j ■ ■ ■ 



x d (l - \x\)H^_ j (W^ 0k _ J ,i J ) C V*(W 7l _0, 



so that the elements of Xd+2-j • ■ 4 %d(^ — \x\)H^_j(W 1 ,o k _ j ,i :i ) satisfy the equation 
(|4.3[) . Secondly, we consider the case of xj+i-j ■ ■ ■ x d H d _.{W~ li o k _ j _ lt x j ,o)- Similar 
to the first CclSC , clS cL result of Corollary 14.31 and Lemma 12.111 we can write 

Xd+i-j ■ ■• x d H^_ j {W lt0k _ j _ uljfi ) C V^(W 7 _iJ, 



so that these polynomials satisfy the equation (|4.3|) for j = 1, 2, . . . , k — 1. Finally, 
the elements of the third term H^(W lt o k ) in the decomposition of W^(W 7 , — lfe) 
satisfy the equation (|4.3[) according to Lemma [2.111 

The same proof applies to the case of k = d + 1 for n > 1 , whereas the case 
Uf (W_i) has been explained right below the statement of Theorem 13.51 

It is remain to prove that the U^iW^ \ k ) has full dimension. From Lemma 
12.101 the restriction of Hf l (W 1 ^ 0k _._ 1 .y) on the face T^~[^~^' is isomorphic to 

Vt^Hw^), so that 

+ d- (k - j + 1) 



dimf^W^o^.,!^) 
which implies, together with dimV^Wy) = that 

fc " 1 'k 



3=1 

AAWu + d-Jt-A/ti + M 
p^Vi/V / V n 

where the last step uses the well-known combinatorial identity. □ 



4.3. Proof of (i) in Theorem 13.71 To show that each bilinear form is an inner 
product, we only need to show that it preserves positivity. Since all coefficients are 
assumed to be non-negative, we clearly have (/, /) _ x > and it remains to show 
that (/, /) _j = implies that / = 0. 

Assume that k — 1 and (/,/)„ _ 1 — 0. Then, by the definition, J^- = for 
1 < i < d, which implies that / is a constant. Furthermore, the last term of 
(/) f) 1 -i shows then / = 0. Thus, (/, /) _ 1 is an inner product. 

Now assume that 2 < k < d and (/, /) _ x = 0. Then each term of in the right 
hand side of (/, /) _j is zero. There are four main terms. The first term shows 

Q x m k = 0, which implies that / has the form 

/(*)= E M x ->)> xj&T d j-\ 

JCm k 
\J\=k-2 
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The second term with | J| = k — 2 shows d dx J — 0, which implies that / has the 
form 

/(*)= E ffc*)' x J^ T Jc 2 ' 

JCm k 
|.7| = fe-3 

we can continue this process with \J\ = j for j = k — 3, . . . , 1 and deduce from the 
second term of (/, /) _ t that / has the form f(x) = f\{x\, . . . ,Xd~k+i)- Now, 
the third term of (/, /} shows that = for j = 1, . . . , d— k + 1, from which 
it follows then that / is a constant. Finally, the fourth term of (/, /) _ x shows 
that f(x) = if k < d or if k = d and A > 0. Thus (/, g) is an inner product. 
The same proof clearly applies to (/, g)_ x - □ 

4.4. Proof of (ii) in Theorem 13.71 We need to show that, for each k, the poly- 
nomial spaces of &<J*(W 7) _ lfc ) given in Theorem 13.51 are the spaces of orthogonal 
polynomials with respect to the inner product (■,■)_, _j . The proof is involved 
and will be divided into several cases. Throughout the proof, we let g be a generic 
element in n~^_ 1 . 

Case 1. k = 1. The inner product (/, g) _ x contains two terms, 

(/,0> 7 ,_i = [f,9]i + [f,9h, 



where 



[/>s]i ■= Y] [ Xi^-^-x 1 dx 

' J T d OXi OXi 



[/,»]»:= A / f(x)g(xW d d _- 1 1 (l-\^i\Y 



i dx < i-i. 

L {d+1} 

By Theorem I3.5[ the space of Z^*(W 7 ,_i) with ji > — 1, 1 < i < d, is decomposed 
into two terms as shown in (|3.4|) . The proof is divided into two cases accordingly. 

Case 1.1. Let f(x) = (1 - \x\)Pi~< A) (x), where P^ A \x) G V^W^i) and n e Nq 
and |n| = n — 1. Then the restriction of / on \x\ = 1 is zero so that [f,g]2 = for 
all g. Furthermore, an integration by parts gives 

(4-5) [f,gh = -flJ Td ~ M)^ 7 ' 1 ^) + (7, + 

Since the term in the curly bracket is a polynomial of degree at most n — 2, the 
orthogonality of Pn 7,1 ^(a;) G V^_ 1 (W 7 ,i) implies immediately that [f,g]i — 0, so 
that (/, g) 1 = in this case. 

Case 1.2. Let / G fl^(W 7)0 ). By LemmaHHUl the restriction of H*(W 1<0 ) on the 
plane \x\ = 1 is exactly V^ _1 (W 7 ). It follows immediately that [/,<?] 2 = for all 
g 6 Hf l _ 1 . Furthermore, an integration by parts gives 

i/.»h = + (■» + 

Since / G V^(W 7 , ) by Lemma[HT0l it follows that [/, g]i = 0, so that (/, g) = 0. 
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By the decomposition of l^(W 7 ,_i), the two cases imply that (/, g) _ 1 = for 
all / £ W„(W 7l _i) and g £ n„_ x , which completes the proof in this case. 

Case 2. 2 < k < d. The inner product (/,<?) _ lfc contains four terms, 
(4.6) (/, g)^_ x = [/, g] x + [/, g} 2 + [/, g]a + if, g]i 

where 

afc— 1 f Pkk — l 



r ak-lfak-l n 



IT 

k-2 



|j|=< 

i=l m k 

[f,9U-=X / /(x d _ fc+ i,0) 5f(x d -fc+i,0)x2_ fc+1 dx d _ fc , 

Jrpd-k 

m k 

where = {d — fc + 2, . . . , d} and = {d — k + 2, . . . , d + 1}. By Theorem l3.51 
the space i/„(W 7) _i fe ), 7, > — 1 for 1 < i < d + 1 — k, is decomposed into three 
pieces as in ()3.4|) . where the second term is a large sum. The proof is divided into 
several cases according to this decomposition. 

Case 2.1. / is an clement of Xd-k+2 4 • ■ Xd+iV^_ k (W J: x k ), the first term of (|3.4I) . 
Let f(x) = • • ■x d x d+1 P^' lk \x), where x d +i = (1 - \x\) and Pn' lk \x) £ 

V*_ fc (W 7lU ) with neNjj and |n| = n - k. 

For I C nifc with |/| = i, the variables on indices in I c = \ / appear as 

product factors in so that the restriction of on the face xfr k+l+1 is zero; 
consequently, [/, g) 2 = 0. For 1 < i < d— k + 1, the variables x d -k+2, ■ • - > %d appear 
as product factors in J^-, so that vanishes on the face T^~ k k+1 ; consequently 
[/iff] 3 = 0- Furthermore, since / vanishes whenever one of x d ~k+2, ■ ■ ■ , x d , x d+ \ 
is zero, it follows that [f,g]i = 0. Thus, it remains to consider [f,g]i- From the 
Rodrigue formula ([21?]) of p£*' lk \x), it follows that 

/(*) = ^-Vif^ {^ U - l)+n (l " kl) |n|+1 } ■ 

Taking derivative of / with respect to variables whose index is in m^, we obtain 

n— 1 f 

- ■ ^y {^>(l-N) W+1 } 

= (i-ki)- (fe - 2) ^i:ur 2)) (-)' 



a*- 1 / _ d k ~ l f _ _ 7 a™- 1 

cb m * ~ dx d - k+ 2 ■••x d _Xd - fc+1 9x n +( c 



from which we obtain immediately that 



f)k-l f f)k-l n 



^ x d+l-fe' r n+(0,l fe _ 1 ) W 9;E m fc 
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Since P^?q^ 2 ^ ( x ) is a polynomial of degree (n — 1) and 

d k ~ 1 q 

n - 1 > (k - 2) + n — k — (k ~ 2) + deg - 

for g £ n^_ l5 it follows from Lemma |2~51 that = 0. This completes the proof 

in this case. 

Case 2.2. / is an element of Xj : k(<j)H^_j (W / 7i ( 0fc _ J .i j ) (T ), where 1 < j < k — 1 and 
c 6 and a^fc contains the variable Xd+i. 

First of all, since the inner product (f,g) rylh involves symmetric derivatives 
with respect to variables {xd-k+2, ■ • ■ , x d }, it is sufficient to consider the case that 
a = id. Let j be fixed, 1 < j ' < k — 1. Since x^fc contains the variable Xd+i, we can 
assume then 

Xj,k = Xd+2-j ■ ■ ■ XdXd+i and /(x) = x d+2 - 3 ■ ■ ■ x d x d+1 Q n ^j(x) 

where Q n -j S H^_j(W 7y o k _ j ^ j ). From Definition 12.91 with Sj being the set Skj — 
{d— fc + 2, . . . , d + 1 — j}, we have then 

f«-i(W r 7.o*-i,i J ) = {^n 7 ' ^'^^) :n£< |n| =n- j withn, =0,i£ 

To keep the notation compact, we introduce the following notations in the proof. 
Recall that for x £ M. d , n £ Nq and 1 < i < d, x$ = (xi,...,Xj) and n; = 
(rii, . . . , rii). For 1 < i < d, we set 

5c j := (x if . . . , x d ) and taj := (rij, . . . , n d ). 

Now, for 1 < j < k — 1, we introduce the notation 

lk,j '■= {l,°k-j,lj-i) and n fcj := (n d -k+x,0k-j,'^d+2-j)- 

Using these notations, we can write the element Q n -j £ S^-j'C^cOfc-iAj) as 

Qn-iCs) - x"^(l - kl)" 1 ^- "• M - |x|)l n «l +1 } , 

where \nk,j \ = n — j. Consequently, our / is of the form 

f)n— j , ^ 

(4.7) f{x) = x d -_\ +1 ^- " M - W)l"^l +1 }. 

The expression shows that /(x) = whenever one of Xd+2-j, ■ ■ ■ } Xd, x d +i is 
zero. It follows readily that [f,g]i = 0, since the indices of {xd+2~j, ■ ■ ■ ,x d ,x d +i} 
are in for 1 < j < k — 1. The variables in {xd+2—j, ■ • • , x d } whose indices are 
in nife appear as product factors in J^, so that vanishes on the face T^~ k k+1 
except when j = 1; consequently [/,<?] 3 = except when j = 1. If j = 1 then 
f(x) — Xd+iQ n -i( x )i where Q n -\(x) £ H^_ 1 (W lt o k _ lt i), so that an integration by 
parts shows that 

d ~ k ~ ir ^ f ( Qg 3^ g 1 

[f,g} 3 = - Al 7 T ™ x 2 + i-^){(7 i + l)^+^^}rfx d - fc+1 . 

i=l m k ^ 

Since / = (1 - M)Q„-i(x) and from LemmafSH Q„_i| T ,-, +1 £ V£lf +1 (W 7 ,i), 

it follows that [/, g]% = in the case of j = 1 as the term in the curly bracket is a 
polynomial of degree at most n — 2. Thus, [/, g\z = for 1 < j < k — 1. 
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Next we consider [/, g] i . For convenience, we introduce the following notation: 
For a fixed j with 1 < j < k — f , define 

m feJ r := {d-k + 2,...,d- j + 1} and m% d := {d - j + 2, . . . , d} 

so that m.k = mfej-Um^ .. In particular, if j = 1, then j = and m£ . = 0. If 

Qfc — 1 

n < k, then |m.fc| = fc — l>ri — 1 and g zm ° = 0, so that [f,g]i =0. If n > k, then 
|m.A;| <n—l. Assume first 2 < j < k—1. Then |m£ .| > 1. We perform integration 
by parts |m£ -| times, once for each of the variables x d +2-j, ■ ■ ■ jX d , which appear 
as product factors in /, to obtain 



(-l)l m S,J 



J Td d+1 ~ k dx" 1 ^ dx m ^ \ dx m " y 1 u J 



Using the product rule, for a generic function h, 

(4.8) ^ {h{x)(l - a;)"} - (1 - a;)"" 1 {(1 - x)h'(x) - nh(x)} , 

it is easy to see that we can write 



where /i n _i mfc i_ i is a polynomial of degree (n — |mfc| — 1). Consequently, we have 

3|m fc . 3 | J 



Now, taking derivative of / in (|4. 7[) with respect to variables whose indices are in 
nifej, we obtain, as these derivatives only apply to the factor (1 — \x\) in /, that 

(4-9) ^— ^ = c K3 x d+2 . 3 ■ ■ ■ x d (l - ^-^-Dp^,!^, -(\ mkJ \-i)) (x ^ 

where Ckj is a constant, x^j = (xd-fc+2) • • ■ ,%d-j+i) and n = n^j. Substituting it 
into the expression for [/, g] 1 , we obtain 

[f>9]i =(- 1 )' m "' jl Cfc,i / xl +1 _ k x d+ 2-j ---Xd 

JT d 

Since |m£ -| = j — 1 and |n| = n — J, we can write |n| — n — .| — 1. As 
|m fc | = |m fcj | + |m^.|, we see that 

A„<r ph'°'-)> 1 i-i-(l m '.il^ 1 )) - (\ m I 1\ 1 A„„h 

deg^n > (,| m fc,il ~ J-J + de g"n-|mj,|-l) 

so that, by Lemma l2~3l [/, g]i = for all g G n^_ 1 when 2 < j < k — 1. In the case 
j = 1, |m|j-| = and there is no need for integration by parts in [f,g]i- In this 
case, (14.91) becomes 

|^ = c M (l - ^D-d^l-Dp^-.-CI^I- 1 ))^) 
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and |n| = n — 1, so that [/, g]i becomes 

jk — 1 £ gk—1 



Since |n| > (|m.fe | — 1) + deg §-g?^ for g G H^^, it follows from Lemma T2.3I that 
[f,g]i = 0. This completes the proof that [f,g]i = 0. 

Next we deal with [f,g]2 — 0. In the proof, let I = I\ U I2, I\ C mfc,j and 
J2 C mjr ■. Firstly, assume that 7 2 is a proper subset of mjr ^ for 2 < j < k — 1. 
Then at least one of the variables in {2^+2-.;, • ■ • ,Xd} whose indices are in m£ . 

appears as a product factor in so that the function ^£ vanishes on the face 
rpd- anc j^ consequently, [/,<7]a = 0. Secondly, assume that /2 = m %,j- Then 
|/2 1 = l m fc j I = j — 1 for 1 < j < k — 1 and the sum over i in [/, 5)2 starts at i = j — 1 
if j > 1, so that we can write 

fc-2 

(4.10) [/, 5 ] 2 = ^ ^ A, /.,/ >.,.,. 

i-raax{l,j-l} 7iCm fc j 

J|=i 

where 

[/>SW : = / T/ „ fc+i+1 00(* - Mp-^-kdxa-k+idxt. 

For j > 1, \h\ = i — 1 > 0. If j > 1, we apply integration by parts IJ2I times 
over the variables Xd+2-j, ■ ■ • ,Xd, which appear as product factors in / by l|4.7j) . to 
obtain 

~ ( ' d+1 - k dx 1 * dx*' \ dx 1 [ lXl) J dXd - fe+lrfX/ 

for i = I J| > j — 1. If /1 t^0 then, using the product rule (|4.8|) . we can write 



■(l-\x\)W- 1 } = (l-\x\)W- 1 h n -\ Ihl , 



dx l2 \ dx 1 

where /i n _m_ 1 is a polynomial of degree (n — |/| — 1), we further deduce that 

[f,g}2, i , I = (-i) |/21 J Td k+i+i *In-*^r {(1 - liD^i-X-m-i} rfx d _ fc+1 d x/ , 

which also holds if j = 1, for which there is no need to take derivatives as I/2I = 0. 
Furthermore, similar to (|4.9[) , if we take derivatives of (|4.7p with respect to variables 
whose indices are in Ii, then the derivatives apply only on the factor (1 — |x|), so 
that 

(4.11) = dujxw-J ■ ■ ■ ^(1 - W )-(IA|-i)^7.o fc - J .i J - 1 ,-(|/ 1 |-i)) (jB)) 
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where n = nfej and dk,j is a constant, which allows us to write 

[f,gh,i,i = (-i) |/2 '4,i / x-d+i-k x <i+2-j ■ ■ ■ x d 



d-k + i + l 



xPi 7 ' 0fc -- 1 - 1 »- (|il| - 1)) (a ; )/ ln _| J |_ 1 ( a; )dx (i _ fe+1( ix z . 

By LemmaElOl the restriction of p^-i' 1 ^ "(mi-i))^ on the face lfr k+i+1 is 
the polynomial Pn 7 ' 0l_J+1,lj_1 ' ^(xd-fc+i, xj) in the (d — + i + l)-variables 

of Xd-fe+i and x/ on 7j c +i+1 . Since |Xa j = j — 1 and |n| = n — j, we have 
|n| > (|7i| — 1) + deg/i and, consequently, by Lemma l2~3l that [/, g]2,i,l — 0. As a 
result, we have shown that the first term of [/, g\i in (|4. 10|) is zero for 1 < j < k — 2 
and 7i ^ 0. 

If, however, Ii = 0, then I = h = m£ ■ , so that i = |/| = j — 1 and, in particular, 
as i > 1, j > 2. We apply integration by parts (|7| — 1) times over the variables 
Xd+3-j, ■■■,%d to obtain 

if &f 
[f,9h,j-i,i = / x 7 fe — /i nHJ |_ l( 2x d _ fc+ idx z , 

where /i„-|/|-i is the polynomial of degree (n — |/| — 1) given by 



0m- 1 ram 



l-MH-1 



dx d+3 -j ■ ■ ■ dx d \ dx 1 
Now, if j > 3, taking derivative of (|4.7[) with respect to a;d+2— j gives 

(4.12) -Of— = X d+ 3-i ■ ■ ■ XdPn 7 ' 0k ~ j + 1 ' lj ~ 2 '°\x), 

dXd+2-j 

where n = n/.^ + ed+2-j and |n| = n — \I\, from which follows immediately that 



[/, g] 2 ,j-i,i = (-l) m 1 J Td _ ik _ j} x] +1 _ fc x d+3 - 



-j %d 

I rpCL— • 

I C 

xPi 7 ' 0fc - +1 ' 1 - 2 ' 0) (x)^_| J |_ 1 (x)rfx d _ fc+1 dx / . 

Since |n| = n — 7 and, by Lemma [2. 101 Pn 7 '° fc ~ 3+1 ' l3 ~ 2 '°' ) (a:)Ld-(fc-j) is an element 
of V^Zm _J (^7,0,1^-3,0) i n the variables (x^-fc+i, xj), it follows from the orthogo- 
nality of V^I^i"^ (W^ai.-^o) that [/)ff]2,i-i,/ = for j > 3. If, however, j = 2, 
then I2 = {d} and [/, g\2,j-i,i for J = 2 becomes 

J1 {d-k+2 ° a 

— df _ p(7,Ojb-i,0), 



By fljjjj) with j = 2, ^ = PA 7 ' J (a;), which, by Lemma is an element 

pd-k+2 

• {d-k+2,...,d-l}' 



in V,^_i +2 (VK 7 ,o,o) when restricted to the face Tf d _^ 2 d _ 1 y Consequently, since 



|n| =?i — l in this case, it follows from the orthogonality of V^Z_\ +2 (W 7) o,o) that 
[/, g]2,i,i = 0. Thus, we have establish that [/, g]2, i, 1 = whenever j — 1 < i < k — 2 
and j > 2. 

Putting these together, we have proved that (/, g} 7 _ lfc = for all g G H^-v 
Consequently, the proof of Case 2.2 is completed. 
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Case 2.3. f is an element of Xj t k(cr)H^_j (W / 7 .( 0fc _ i .ij)cr) , where 1 < j < k — 1 and 
c 6 and Xj^fc does not contain the variable Xd+i- 

As in the Case 2.2, it is sufficient to consider a = id. Since Xj k does not contain 
Xd+i, we can assume that 

Xj,k = Xd+2- 3 ■■■Xd and f(x) = x d+ 2- 3 ■ ■ ■ x d Q n -j+i(x) 

for 2 < j < k, where Q n -j+i £ H^_ j+1 {W 1<0k _ j! x j _ 1 fl). From Definition El 
x s k - j+1 = {1 - \x\,Xd-k+2, ■ ■ ■ , Xd-j+i} and then 

H^W y , Oh _ iflj _ u0 ) = {pi 0k - i+1 ' 1 ^\x Sk _ j+1 ,*d- j+ 2,Xd-k) : 

n = (0, n fcj ), |n| = n, n k j £ N^^^ 1 } . 

From Definition 12.41 we can explicitly write the Rodrigue basis of this space. To 
keep the notation compact, we use following notations in this case of the proof, 

x fc ,j := (£ d - j+2 ,Xd-k) and -y kJ := (lj-i,7d-k)- 
In case 2 of Definition 12.41 taking 

Xi 2 X d — k-\-2-> ■ • ■ ; Xi k — Xd, Xi k ^ 1 — Xi, . . . , Xi d+1 — Xd—k+1 

n = (n k -j+i, n kij ) with n k - j+ i = (m, . . . , nu-j+i), n k j e N ~ 
then we can write 



- f-l^iY^J'T -7 ''-^ 1 -- f„ n '=-j + i Y "''fc.j+ n fc.i T ,7d~fc+l + |nn 

where diikj = ^d-fc+2,d-fe+i i ■ • • > d d ,d-k+i, di,d-k+\, • ■ • , ^d-fc,d-fc+i) and 

recall that 9jj = di — dj. Since the first {k — j + 1) components of n are zero, in the 
definition of H^(W 1 ^ k _^ 1] _ 1 ^), we can write Q n -j+i & H^ j+1 (W liOk _ jAj _ u0 ) as 

n ■ , ,M - v" 7t -' T " 7i - |,+1 ^ |nfc ' jl / v ' Y *' ;,+n *' i T 7 ' i -* +1+|,1 *^ |, l 

Vn-J + lW- x );j A <i-fc+l a n k ,j \^k,j x d-k+l ] > 

where <9y k j = (9<i-j+2,d-fe+i, • • • , 9d,d-fc+i, <9i,d-fc+i> . . . , d d -k,d-k+i)- Hence, we 
conclude that 



(A f^-l - „-Td-fc + i ^" 3 + 1 f Tfc.j+^.j |nfc, J l+7d-fc + ll 

^• i0 ^ J W — x d-fc+l Qy n k,j \ k,j x d-k+l j ■ 



where |nfc,j| = n — j + 1. 

This expression of / shows that f(x) = whenever one of Xd+2-j, ■ ■ ■ , Xd is zero. 
It follows readily that [f,g]i = 0, since the indices of {x d +2—jj ■ ■ ■ ,Xd} are among 
the variables whose indices are in for 2 < j < k, so that the restriction of / 
on the face T d + fc is zero. Furthermore, [f,g]s — for the same reason. Since the 

variables in {x d -k+2, ■ ■ • , x d -j+i}, 2 < j < k - 1, do not appear in /, = 0, 

so that [/, g\i — except for j — k. We now consider the case of j = k. In this 
case, by (|4.13|) with j — k, 

an—k+l . I,, 
(4 141 f(r) - x^-fc+i / x Tfc.*+ n *.fc Jn*,fc|+7«(-* + il 

/ — X d-k+l Q y n k,k \ X k,k X d-k+l \ ■ 



k 



SOBOLEV ORTHOGONAL POLYNOMIALS ON A SIMPLEX 



29 



where | , fc | = n — k + 1. Applying integration by parts |mfc| — 1 times over the 
variables Xd-k+3, • ■ • , Xd, which are product factors in /, we obtain 

[/,.,], - (-1)'-"- L*u-* JlJZl - w)lm ""'} d * 

where M — = - — - — H — . Taking derivative of f in (14.141) with respect to 

Xd-k+2, we obtain 

3f D(0,0,lfe_2,7)/ \ 

- CX d -k+3 ■ ■ -XdPa y (^S fe _ 3+1 ,Xfe,j), 



dxd-k+2 

where c is a constant and n = (Q, n*,,*,). Since g Xd hl k+3 { g^mf (1 — |x|)' mfc ' _1 | is 
a polynomial of degree n — |mfc| — 1 and |n| — n — k + 1 = n — |ixifc|, we deduce 
that = form the fact that p^ ' 1 *-*^ g V^_ |mfc| (W (7 , ,i fe _ 2 ,o))- Thus, 

[/, 9 ] 1 =0for2< j<fc. 

It remains to consider [/, g\2- Firstly, assume that I2 is a proper subset of m£ ■ 
for 2 < j < k. Since at least one of the variables in {xd+%— j, • ■ • , Xd} whose indices 
are in m c k j appears as a product factor in , the function vanishes on the face 
T d-k+i+x g0 that ^ g j 2 = q Secondly, assume that I 2 = m L k j = {d- 3 + 2, . . . , d}. 
Then | 1 = |ni|j| = j — 1 for 2 < 3 < k — 1. We do not need to consider the case 
j = k when I2 = ., since if j = k then m^j = and = m c k ■ = I2, so that 
/ = rrifc, which contradicts to the fact that / is a proper subset of in [f,g}2- 
Thus, we only need to consider 2 < j < k — 1, for which the summation over i in 
[/, g]2 starts at i = 3 — 1, so that we can write 

fe-2 

(4-15) [/, ff ] 2 = £ ^ A 7 [/, 5 ] 2 , a , 

i=j-l IiCmtj 
|J|=i 

where [/, 5)2,1,1 are given as in (14. 10[) . Since the variables Xd-k+2, ■ • - > Xd-j+i whose 
indices are in m^- for 2 < j < k—1 do not appear in /, dx i{ — and, consequently, 
[/, ff] 2 ,i,i = if i > j. For i = j-1, \I 2 \ = \I\ =3-1 and % = 1= {d-j+2, . . . , d}; 
in particular, /| = / c = {d — k + 2, . . . , d — 3 + 1}- Applying integration by parts 
IJ2I — 1 times over the variables Xd+3-j, ■ ■ • , Xd, which appeared as product factors 
in /, we find 

f df 

[f, 9)2,3-2,1 = / 1 _ fe - — ■ h n -j{x) dxd-k+idy-i, 

JTjc OXd+2-j 

where h n -j is a polynomial of degree n — j given by 

h (x) U ' 1 P' "A 
Now, taking derivative of / in (|4.13|) with respect to Xd+2-j, we obtain 



Of 



1 7~>(0fc — 1 4- 2 ) 1 7 — 2 j T ) / \ 

= cUd-i+a • • -XdPn J+ (a;s fc _ j+1 ,x fcj ), 



dXd+2-j 

: a: 

is, by Lemma [2.10[ an element of V^Zj+'i 1 (W 7 ,o,i^_ 2 ,o) when restricted to the face 



where d is a constant and n = (0, n^j). The polynomial Pn° fc 3+2 ' lj 2,7 ' ) (a;Sfc_j + i 7 x fc,j) 
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d-j+i}- Consequently, since |n| = n — j + 1 and h is a polynomial of de- 
gree n-j, [f,g]2,j-i,i = follows from the orthogonality of V^Zf+i (W 7 , ,i 3 _ 2) o)- 
Hence, [/,5]2 = for all 5 £ n£_x- 

Putting these together, we have proved that (/,</) _ lk = for all g G n^_ 1 . 
Consequently, the proof of Case 2.3 is completed. 

Case 2.1 f G fl^(W 7>0 J. From Definitional we can assume 

f( x ) = Pn'' l) ( x d-k+2, ■ ■ ■ ,x d+ i,x d _ fe ), |n| = n, 

where x d+ i = 1 — \x\ and n = (0, n d _ fc ), n d _ fc = (m, . . . , ra d -fc). Since the first 
k components of n are zero, the Rodrigue formula (|2.9[) shows that the variables 

qk — 1 f Ql f 

whose indices are in do not really appear in /. Hence, dx ™l = and ^p- = 
for I C mfc, from which it follows that [/, 5)1 = and [/, 5)2 = 0, respectively. For 
the integral in [/, 5)3, we apply integration by parts to obtain 

d k- 1 « 1" ^ <9 2 5 1 

J— 1 m fc v J 

Since, by Lemma l2.101 f\ T d-k+i G Vn~ fc+1 (W 7i o), an d *he term in the curly bracket 
is a polynomial of degree at most n — 2, we conclude that [/, 5)3 = 0. Finally, by 
Lemma T2.101 f\ T d-h G V^ _fe (W 7 ), which implies immediately that [/, 5)4 = 0. 

Thus, it follows from gU) that (f,g) _ lfc = for all 5 G n^_ x . 

Putting the four cases together, the decomposition (|3.4[) shows that {/, 5) _ x = 
for every / G ^(W-y.-ijJ an d all 5 G n^_ 1 . This completes the proof of (ii) of 
the theorem for Case 2. 

Case 3. k = d+ 1. The main terms of the inner product (/, g)_ 1 are essentially the 
case of k = d + 1 of [/, 5)1 and [/, 5)2 in Case 2, hence, as shown in that case, that 
these two terms are zero for every / G U„(W-±) and for all 5 G n„_ x . Furthermore, 
if / G x\ ■ ■ ■ x c i+iV^_ d _ 1 (H / i) then / vanishes on eo, ei, . . . , e<j, whereas if / G 
x Jiti +i(o-)^_ J (W / ( o d+1 _ 3 i1)ct ), then /(e ) = /(ei) = ... = f(e d ) = 0. Thus, every 
/ G U^(W-i) vanishes at all vertices, so that the last term in {f,g)_ 1 is zero. 
Consequently, (f,g)_ ± = for V/ G and V5 G IL*_ V 

The proof of (ii) in Theorem 13. 71 is completed. □ 

4.5. Proof of (iii) in Theorem 13. Tl We need to show that (f,g)„ _ 1 — for 
/ G (l-|£|)Kf-i(W 7 ,i) and 5 G i^(W 7 , ). The proof uses a property of i^(W 7 ,o), 
which is of interesting in itself. 

By definition, H%(W 1>0 ) = H% {d+1} (W 7t a), that is, Si = {d+1} in the item 2 of 
Definition and (js,Ji d ) = {0,l',7e d ), where 7' = 7s\{d+i} = ("fd, ■ ■ • 
and a; s = (xSi,xs\Si) = (%d+i,x'), x 1 = (x ei , . .., xe d _ 1 ). Hence, by Definition H31 
the elements of if^(W 7i o) are given by, since ri\ = 0, 

P„ (x s ) = (x) ix ld d —, 

where n = (0,n') and 9xs = (di lt t d , . . . , d£ d _ lt £ d ). Consequently, the elements 
of ff^(W 7: o) are homogeneous polynomials of degree |n|, which satisfy, by Euler's 



) (^)=(,o-v:^{(-'r' +n '^ +|n| }, 
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formula, the equation x\d\g + . . . + Xdddg — ng. Taking derivative again shows 
easily that, if g is a homogeneous polynomial of degree n, then 

d 

^2xfdfg + 2 ^2 XiXjdidjg = n(n-l)g, 

i=l l<i<j<d 

which applies, in particular, to g £ H„(W^ t o). On the other hand, by the item 
(i) of Lemma \2 . 1 01 and Lemma \2. Ill g £ H^(W 7 .o) satisfies the equation Ly,-ig = 
—n(n + |7| + d - l)g, which is, by (fl~2j) 

d 

^ Xi(l - Xi)d\g -2 ^2 XiXjdidjg 

i=l \<i<j<d 
d 

+ ^ + 1 ~ d 7 ' + c ar *) d *9 = -n{n+\j\+d- 1) g. 

i=l 

Adding the two identities together and making use of the Euler's formula again, we 
conclude that g £ H%(Wy,o) satisfies the following equation 

d d 

(4.16) x ^9 + 5^(7i + 1)^.9 = 0. 

1=1 8=1 

Now, let f(x) = (1 - \x\)P^' X) (x), where P^'^ix) £ V^W^i), and g £ 
H*(W~t o)- Following the proof of Case 1 of the subsection 4.4, we then arrive at 
(|4.5p . The summation over the derivatives of g, however, is exactly the differential 
operator in the left hand side of (|4.16p . so that (f,g) _ 1 =Q. □ 
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